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This work is concerned with the instability and 
unbalance response of dissymmetric rotor-bearing systems 
containing periodic coefficients when modeling produces 
matrices with a large number of degrees of freedom. It 
is important to solve the equations and then to predict 
the dynamic behavior of the system. This can be done 
knowing the instability areas, and the unbalance 
response in the stable areas. One deals here with a 
large number of equations and a reduction of the number 
of degrees of freedom of the system is achieved through 
a pseudo modal method. This method is shown to give 
satisfactory results. 

INTRODUCTION 

It is now more and more necessary to predict accurately the dynamic 
behavior of rotor bearings systems : natural frequencies as a function of 
the speed of rotation vAiich gives the critical speeds and instability zones 
and response to unbalance and nonsynchronous forces. 

For most of rotors containing a large number of degrees of freedom, 
calculations do not pose any problem today. Many authors, as £lj , £ 2 ], [V] , 

, have shown that, using a method based either on substructuring or on a 
modal reduction, the number of degrees of freedom can be lowered signifi- 
cately and the results are without significant loss of accuracy. 

For dissymmetric rotors, it happens sometimes that the equations lead 
to periodic coefficients. In this case, some authors, as [ 5 ] , [ 6 ] , [[ 7 ]], 
give methods to detect instability zones (mainly transfer-matrix or 
infinite determinant methods) and to calculate step by step the unbalance 
response. But usually, these methods are applied on systems having a few 
degrees of freedom (up to 20). The purpose of this study is to develop for 
large periodic systems, a method \diich can be used to obtain the dynamical 
behavior of the systems. In a first part, a pseudo-modal method adapted to 
the resolution of periodic differential equations is presented. In a second 
part an industrial application containing 96 degrees of freedom is con- 
sidered. The influence of the number of modes on the results, and the 
accuracy of the methods are discussed. 

PSEUDO-MODAL METHOD 

The systerae to be solved can be written as : 

MX* 0 + A(t) X° + K(t)X = F(t) (1) 
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M, A, K are n order matrices and it has been shown [JT] that A(t) and K(t) 
can be written as : 


A(t) = A + A. sin 2Qt + A 0 cos 2£2t 

o 1 z 

K(t) = K + K, sin 2Qt + L cos 2Slt 

o 1 Z 


( 2 ) 

(3) 


The pseudo-modal method consists in using a "modal" base of the system 

M X* ° + K X = 0 (4) 

The modal base is calculated from the system at rest (Q = 0) and at 
the initial instant (t = 0). In these conditions, there are neither 
periodic coefficients, nor gyroscopic effects. Furthemore , the damping of 
the bearings is omited , and the non symmetric terms of the stiffness matrix 
K are symmetrised : 


K = K q+ K 2 
K = sym (K) 


(5) 

( 6 ) 


The modal base 0 is built with the lowest 2, modes of (4), (2. « n) . 

The relation between the degrees of freedom of the system and the modal 
parameters is : 

X = 0 q (7) 

Equation (1) with (2), (3) and (7) leads to : 

M 0 q** + (A + A, sin 2Qt + A sin 2fit) 0 q* + (K + K sin 2fit 

r ° 1 2 ° 1 ( 8 ) 
+ K 2 cos 2Qt) 0 q = F 

and premultiplying by ft , (8) can be written as : 

m q** + (a +a, sin 2ftt+a„ cos 2J2t)q* + (k +k. sin 2nt+k„ cos 2Qt)q = f (9) 
^ o 1 2 o 1 L 

with 

m = f M 0 ; a = ft A 0 ; a = ft A 0 ; a = ft A 0 

° o l l z z (10) 

k Q = ft K q 0 ; kj = 0 ^0 •, k 2 = 0 K 2 0 ; f = 0 F 

SOLUTION OF THE EQUATIONS 

Instabilities and unbalance response come from the reduced system (9), 
and from (7). The method used is that detailed in []8^ • Here the basic 
principles are shown. 

* Instabilities are found with the resolution of : 


m q* # + a q° + kq = 0 with : 


a = a + a, sin 2Qt + a„ cos 2£2t 

o 1 2 

k = k + k. sin 2Qt + k n cos 2Qt 

o 1 Z 


( 11 ) 

( 12 ) 

(13) 


(11) is transformed into : 
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dp 

= B p 

dt 


with 


P = 


and 


q 

e 

q 


B 



k 



(14) 


(15) 


(16) 


B is a periodic matrix, of period T. The period T of the periodic coeffi- 
cients is divided in s intervals of length h = T/s. System (14) is consi- 
dered to be constant on each interval and a matrix T. connecting displa- 
cements and velocities at the instants ih and (i-l)h caA be calculated. The 
general matrix connecting q and q* over a period T is obtained by the 
product of matrices : 


T, = T - 
f s — 1 


T. . T. . 
J J-l 


. . . T, . T 
1 o 


(17) 


and 


p(T) = T f p(0) 


(18) 


The 2.1 complex eigenvalues of T^ are representative of the stability 
of the system. If they are all less than unity the system is stable. Here 
each matrix T. is calculated with a Newmark formulation and the expression 
of T. is given-*below. 


T. 

J 


with 


D -1 F 


I — (D -1 F-I) 


d" 1 E 


-1 

— D E- I 


(19) 


D = 


E = 


F = 


4m 

4m 

h 

4m 


2c . . 
J + l 


+ k 


j+i 


+ c . . - c . 
J+l J 


2c. , 

J + 1 


- k. 
J 


( 20 ) 


( 21 ) 


( 22 ) 


Instability zones are obviously the same as those which would be 
obtained from (1). 


* Unbalance response is the solution of (9). Numerical resolution is 
made with a Newnark formulation and the initial conditions are chosen as : 

t = 0 q(0) = 0 (23) 
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q* (0) 
q* * (O) 


0 

m * f(0) 

APPLICATION 




(24) 

(25) 


A rotor of a 220 KW Steam Compressor, (Fig.l) is studied. The rotor 
contains symmetric shaft and disks, and dissymmetric bearings. Calculations 
are at first made in a fixed reference frame where equations have constant 
coefficients and are obtained with a finite element model and are easily 
solved £ 93 . Then, calculations are made in a rotating reference frame : the 
equations have periodic coefficients and are also obtained with a finite 
element model. The comparison of the results in the two different frames 
shows the interest of the method proposed. 



Fig. 1 : 220 KW Steam Compressor Rotor 

1 - Description of the model 

The rotor is modeled with 23 finite elements, as shown in Fig. 2. Two 
calculations are made with different materials for the disks. They are in 
Aluminium for rotor A, and in steel for rotor B. 



BEARINGS 



Fig. 2 : Modelisation of the rotor 

2 - Calculations 

The two rotors are used for different purposes. For rotor A (aluminium 
disks) the unbalance response can be performed over 4 critical speeds 
because the motion is always stable in the operating range (0 - 60000 RFM) . 
For rotor B, instabilities appear near 20.000 RPM. So the two configura- 
tions will provide satisfactory tests both for instabilities and unbalance 
response . 


3 - Results 


* Instabilities 

They appear at 21300 RPM in the fixed reference frame. In the rotating 
one the instability with 14 modes appears at 19700 RPM. 

* Unbalance response 

Fig. 3, 4, 5 represent the maximum of the unbalance response at node 
1. Different number of modes (6-10-14) are considered, and compared to the 
results in the fixed-reference frame. The resolution in the rotating 
reference frame with reduced coordinates introduces a slight gap in the 
frequencies. This creates differences in the amplitude to up to 10 % in the 
critical frequencies zones. Only the permanent solution is compared, and 
the number of rotations of the rotor necessary to obtain the permanent 
solution depends a lot on its rotating speed. 

CONCLUSION 

The dynamic behavior of periodic coefficient large system is predicted 
here with a pseudo modal method, using a significant reduction of the 
number of degrees of freedom. 

Differences which are shown in the industrial example presented may be 
inherent to the numerical calculations. 
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NOMENCLATURE 

A "Damping" matrix [n.nl containing periodic coefficients 

a modal damping matrix 

B matrix which transforms a [jl.S,] second order differential 

system into in first order differential system 

E Young modulus (N/m ) 

F unbalance force vector 

h time interval h = t/s 

I unity matrix 

K "stiffness" matrix (jn.nl containing periodic coefficients 

k modal stiffness matrix jjL.jT] 

£ number of modes taken into acount 

M mass matrix [[n.n] 

m modal mass matrix 

n number of degrees of freedom 

P vector containing modal displacements and vitess 

q modal displacement vector ( 

s number of intervals in a period T 

T period of coefficients of differential equations 
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b 

a 


P 

V 


• • 


transfer matrix over one period T 

displacements vector QT| 

modal base built with the first l modes 

rotation speed of ^he rotor 

volumic mass (Kg/m ) 

Poisson coefficient 
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